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The penetration of an ac magnetic signal into a type II superconductor residing in the Shubnikov 
phase depends on the pinning properties of Abrikosov vortices. Within a phenomenological theory, 
the so-called Campbell penetration depth Ac is determined by the curvature a at the bottom of the 
effective pinning potential. Preparing the sample into a critical state, this curvature vanishes and 
the Campbell length formally diverges. We make use of the microscopic expression for the pinning 
force density derived within strong pinning theory and show how flux penetration on top of a critical 
state proceeds in a regular way. 
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I. INTRODUCTION 

Transport 1 and magnetization 2 measurements are well 
known tools for the basic phenomenological characteri¬ 
zation of superconductors. The diamagnetic screening 
in a superconductor involves multiple aspects: a weak 
magnetic field penetrates to the material over the length 
A l , the (London) penetration depth' 3 , which provides ac¬ 
cess to the superfluid density p s and is typically of sub¬ 
micrometer size. In type II superconductors, a magnetic 
field H penetrates the material through quantized flux¬ 
lines or Abrikosov vortices 4 which arrange in a triangular 
lattice defining the Shubnikov phase 5 . Testing this phase 
via a small ac-magnetic field produces a normal response 
described by the skin-effect, with a reduced (flux-flow) 
resistivity ps entering the usual expression for the skin- 
depth 6(to). In real materials, vortices get pinned by 
material defects, thereby establishing the desired criti¬ 
cal current density j c below which vortices are trapped. 
In this situation, an external ac-field probes the pinning 
landscape (or pinscape), as the latter now determines the 
field penetration over the scale A c (typically 1 100 mi¬ 

crometers), the Campbell length 5 . The linear response in 
the Campbell regime assumes that the vortex displace¬ 
ments induced by the ac-magnetic field h ac are much 
smaller than the characteristic pinning length such that 
the Campbell response probes the pinning wells. This 
condition implies that the currents Sj ~ ch ac /\ c induced 
by the vortex displacements are much smaller than the 
system’s critical current j c , Sj <C j c • Increasing the field 
strength h ac to large values capable of changing the di¬ 
rection of the critical state periodically, i.e., Sj > j c , the 
response is described by the Bean model'. In this limit, 
the Bean penetration depth i B ~ ch ac /4-Kj c depends on 
the field amplitude h ac and hence generates a higher har¬ 
monic signal in the magnetic response. 

Within a phenomenological model of the Campbell re¬ 
sponse, the dynamics of vortices in a random potential 
landscape is reduced to the motion of a vortex in an effec¬ 
tive defect potential. Vortices probe the bottom of the 
pinning potential and the Campbell penetration depth 
A c oc 1/i/a involves the curvature a at the minimum 
of the effective pinning well 5 . This description becomes 



FIG. 1. Sketch of the linear Campbell (left panel) vs. non¬ 
linear Bean (right panel) response of a vortex critical state 
subject to an ac magnetic field h ac . The gray-shaded areas 
show the change of the dc field deep inside the sample after 
the initial field penetration at short times. The blue regions 
indicate the magnetic field oscillations at large times pene¬ 
trating the sample to the depths A c and t B , respectively. The 
cross-over between the Campbell- and Bean regimes with in¬ 
creasing field amplitude h ac occurs when the induced current 
5j near the sample surface is of the order of j c , i.e., when 
h ac > jc A c /c. 


quite problematic when dealing with the technologically 
most relevant vortex configuration, the critical state' re¬ 
alizing the maximally possible current flow (the criti¬ 
cal current density j c ) before depinning. Indeed, within 
this approach, the critical state is characterized by a 
vanishing 8 curvature a(j) oc ( j c — j) 1 ^ 2 —> 0 and thus 
implies a diverging Campbell length A c . This result is 
not satisfactory in two respects: first, the predicted di¬ 
vergence is not observed in experiments 8 ; an explanation 
that flux creep prevents probing of the proper critical 
state is hardly applicable to the early work by Campbell 5 
on low-T c Pb-Bi alloys. Second, it turns out that the 
Campbell length in the critical state can be even smaller 
than that in the field-cooled state 9,10 . 

In this paper, we make use of a microscopic the¬ 
ory in order to reconcile the apparent divergence of 
the penetration depth A c oc 1 /yfa in the critical state 
with a regular vortex dynamics. To this end, we de¬ 
scribe the ac-magnetic penetration within the strong pin¬ 
ning framework 10 and determine the dynamical evolution 
of the vortex state when probing a critical state with 
a small-amplitude and low-frequency ac-magnetic field 
h ac e~ iult - We show that in this situation the ac-response 
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involves a transient region where vortices first penetrate 
throughout the sample in the form of diffusive flux pulses. 
The penetrated flux produces an upward shift of the crit¬ 
ical state profile; once this dc shift reaches the maximal 
field amplitude h ac , see Fig. 1 (left), the ac-field only low¬ 
ers the magnetic held at the sample edge and we arrive at 
a finite ac penetration depth A c oc l/-^/A/ p i n , with A/ pin 
the jump in the microscopic pinning force for the criti¬ 
cal state as given by the strong pinning theory. While, 
qualitatively, a similar picture describes the non-linear 
ac response at large held amplitudes h ac > j c X c /c as de¬ 
scribed by the Bean model 7 and illustrated in the right 
panel of Fig. 1, the quantitative description and results 
are very different for the Campbell and Bean penetration 
regimes. 

In the following, we briehy review (Sec. II) the ac- 
magnetic response of a type II superconductor in the 
Shubnikov phase and its microscopic extension using the 
result from strong pinning theory. In Section III we de¬ 
scribe the transient regime of the ac-response with its 
vortex penetration, the central topic of this paper. Sec¬ 
tion IV gives a short summary and conclusions. 


II. ac RESPONSE AND CAMPBELL LENGTH 

An ac dynamical held enters a homogeneous metallic 
sample over the skin-depth S(uj) ~ \Jc 2 /2nfj,u>a, with 
/i and a the materials’ magnetic permeability and con¬ 
ductivity, respectively. For the inhomogeneous state of a 
held-penetrated type II superconductor, it is the behavior 
of vortices which determines the ac response. In a real, 
i.e., defected material, the dynamics of pinned vortices 
is dictated by the pinning potential opposing their mo¬ 
tion; within a phenomenological model, the force density 
F pm ~ — ax dehnes the penetration depth A c at small 
frequencies u>, 


A 


2 _ 
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where B is the mean penetrated held. The above re¬ 
sult for A c was hrst given by Archie Campbell 6 , together 
with supporting experimental data on the ac-magnetic 
response of a superconductor. 

For a brief derivation of the result (1), we consider a su¬ 
perconductor occupying the half-space X > 0 in the pres¬ 
ence of a magnetic held involving both dc and (a small) ac 
component H(t) = Hq + h ac exp(— iwt) directed along Z 
(for consistency with Ref. 10 we use capital-letter coordi¬ 
nates in describing the macroscopic situation). This held 
penetrates to the sample in the form of vortices producing 
an average magnetic induction B ( X , t). The screening 
current density j in the superconductor hows along the 
E-axis and exerts a Lorentz force density F L = jB/c di¬ 
rected into the sample. In a stationary state, the Lorentz 
force density F, has to be balanced by the pinning force 
density F pin , otherwise vortices move dissipatively. In¬ 
troducing the macroscopic displacement held U(X,t) of 


the vortex system, the force balance equation takes the 
form 


rjd t U = F4j,U) + F pin (X,U), (2) 

with p denoting the viscosity 11 . The magnetic induction 
and current can be split into a dc part and a contribution 
from the external ac drive, B(X,t) = B 0 + SB(X,t) and 
j(X , t) = jo+Sj(X, f), where B is driven at the boundary, 
B(0,t) = Bo+h ac exp(—iLjt). The ac-magnetic induction 
SB relates to Sj via Ampere’s law, d x SB = — (c/47t )Sj, 
and to the displacement U via the change in vortex den¬ 
sity, SB/Bq « —d x U. For a critical state, the dc current 
density jo is maximal and hence equals the critical cur¬ 
rent density j c . The latter is compensated by the max¬ 
imal (or critical) pinning force density maxF p j n = — F c , 
with F c = j c B/c. Rewriting the Lorentz force density 
through U and denoting deviations from maximal pin¬ 
ning by SF pm = F c — F p i n (U), we arrive at the dynamical 
equation of the form, 

V&tU - (B 2 /4Tr)d 2 x U - 6F pin (U) = 0, (3) 

Alternatively, this equation can be obtained starting from 
the equation of motion of individual vortices. Aver¬ 
aging over many inter-vortex spacings ao, one then ar¬ 
rives at the above expression. With the chosen field and 
current directions, the vortex-vortex interaction [second 
term in Eq. (3)] only involves the bulk compression mod¬ 
ulus cn(fe = 0) = Rq/47t, while contributions of the 
shear and tilt moduli are averaged to zero. 

The non-trivial part in arriving at an explicit equa¬ 
tion for U(X, t) is the functional dependence of the pin¬ 
ning force density SF p i n (U). Assuming vortices trapped 
in a harmonic pinning potential, Campbell 6 introduced 
the phenomenological Ansatz SF pin (U) = —aU , with a 
describing the curvature of the effective pinning poten¬ 
tial. The resulting differential equation is of the driven- 
diffusive type and easily solved, 

U(X, t) = \ c {h ac /Bo)e~ x ^ Xc e~ iut (4) 


with 


A» 


B% 1 

4-7T a — icorj 


( 5 ) 


At small frequencies, we obtain the Campbell length 
Ac = A c (u> = 0) as given in Eq. (1); at high frequen¬ 
cies we make use of the Bardeen formula r/ = BH&fc 2 p n 
to arrive at the skin depth 8(u) ~ (c 2 /27ra;crff) 1 / 2 with 
ag = cr„H C 2 /B the flux flow conductivity, a n the nor¬ 
mal state conductivity, H c2 = $ 0 /27t£ 2 the upper critical 
field, and £ the coherence length of the superconductor. 

Here, we go beyond this phenomenological theory 
and make use of the expression for the restoring force 
SF pm derived from strong pinning theory 10 . The lat¬ 
ter goes back to early work of Labusch 12 and of Larkin 
and Ovchinnikov 13 and has attracted quite some inter¬ 
est over the recent years 14-16 . Within the framework 
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FIG. 2. Microscopic effective pinning force / p i n with bistable 
pinned and unpinned solutions as a function of the pin-vortex 
distance x (dotted lines denote unstable solutions). For ran¬ 
domly positioned defects and the vortex lattice in the critical 
state, the occupation of the force branches (blue) produces 
the maximal restoring force F pln = n p (f pln ) = —F c . A macro¬ 
scopic displacement U > 0 of all vortices in the direction of 
the Lorentz force does not change the branch occupation and 
8F p i n = 0. A displacement U < 0 against the Lorentz force, 
however, modifies the branch occupation (green) and reduces 
the pinning force SF pin oc Af pin U, see Eq. (7). 


of strong pinning theory, vortices are pinned by indi¬ 
vidual defects, thus allowing for a quantitative descrip¬ 
tion of pinning related phenomena. The crucial feature 
appearing within strong pinning is the strong deforma¬ 
tion of vortices giving rise to bistable solutions of the 
force equation balancing the elastic vortex energy against 
the pinning energy due to the defect; the appearance 
of such bistable solutions, quantitatively formulated in 
the Labusch criterion 12 k = 1, then separates strong 
(k > 1) from weak (k < 1) pinning (the Labusch pa¬ 
rameter n ~ fp/Ct; measures the relative strength of the 
pinning force f p of one defect as compared to an effective 
elasticity C of the vortex lattice). 

Within strong pinning theory, one studies how a rep¬ 
resentative vortex embedded in the vortex lattice gets 
locally deformed and pinned due to the presence of a 
defect 10 ; the macroscopic pinning force density F p i n then 
results from proper averaging of the microscopic effec¬ 
tive pinning force f pin . More specifically, when dragging 
a vortex across an individual defect, the vortex jumps 
into the pinning potential at — x _, X- ~ £, and remains 
pinned therein until the deformation becomes too large 
and the vortex snaps out of the pin at x+ ~ k£, see 
Fig. 2. Assuming a small density n p of pinning centers, 
one can ignore interactions between defects and the pin¬ 
ning force density F p j n = n p (f p i n ) derives from a simple 
average over the microscopic pinning states with vortices 
occupying pinned and unpinned branches of the bistable 
pinning landscape. A vortex system in the (Bean) criti¬ 
cal state is defined through the maximal averaged pinning 
force density F c produced by the pins. In this state, the 
vortex configuration is critical, i.e., when shifting the vor¬ 
tex system in the direction of vortex penetration U > 0 


there is no change in pinning force as the latter is already 
maximal, 

SF pin (U > 0) = 0. (6) 

On the other hand, moving vortices opposite to the crit¬ 
ical slope, i.e., for U < 0, the branch occupation rear¬ 
ranges, see Fig. 2, and the pinning force density is di¬ 
minished by 10 

6F pin (U < 0) » -n p (t ± /a 2 0 )Af pin U. (7) 

Here, A/ p ; n denotes the sum of jumps in the pinning 
force when vortices jump into and snap out of the pin¬ 
ning trap created by a defect. Furthermore, t± is the 
transverse length over which vortices passing by the de¬ 
fect are trapped and a p 2 = -B/$o is the vortex density 
(<f>o = he/ 2e denotes the flux unit). For a point-like de¬ 
fect, the transverse trapping length is of the order of the 
vortex core size, t± ~ £. 

The drop (7) in critical force appears whenever vortices 
start moving to the left, i.e., when U(X, t) decreases with 
increasing t. In the dynamical situation defined by the 
ac-response, we have to follow the macroscopic displace¬ 
ment U ( X , t) in time and switch on the restoring force 
(7) when U starts decreasing. Assume that vortices have 
reached the displacement Uq when they change direction 
of motion; the argument in Eq. (7) then has to be re¬ 
placed by U—Uo, with the restoring force smoothly grow¬ 
ing from zero. In a fully dynamical situation, Uq has to be 
calculated self-consistently and is given by the maximal 
displacement reached so far, Uq(X, t) = max t / <t U(X, t'). 
The final expression for the reduction in pinning force 
then is given by 

6F pin (U) = -a sp (U-U 0 ) ( 8 ) 

with 

C^sp ~ ( 9 ) 

In the following, we solve the dynamical equation for 
U(X,t ), Eq. (3), with the restoring force SF p in , Eq. (8), 
derived from strong pinning theory. 

III. CRITICAL STATE ac-RESPONSE 

Before entering the detailed discussion we give a short 
overview on the ac dynamics. Inserting the result (8) into 
Eq. (3) generates a complex vortex dynamics as flux en¬ 
ters the sample in a sequence of diffusive pulses until the 
internal magnetic field is raised to B 0 + h ac —the discus¬ 
sion of this initial dynamics is the central topic of the pa¬ 
per carried out below. Once the asymptotic time domain 
has been reached, vortices exhibit the typical oscillatory 
behavior within the pinning wells, but with respect to 
the new critical state that has been shifted upward by 
h ac , see Sec. IIIB below. The ac magnetic response then 
follows the standard result with the Campbell length A c 
determined by the jump in pinning force A/ pin through 
a sp , see Eq. (9). 
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FIG. 3. Displacement U (solid, blue) and max-field Uo (black, 
dashed) as a function of time t at the positions X = 0 (sam¬ 
ple boundary), X = 2 Ac, and X = 4Ac (subsequent curves 
are shifted for better visibility). The periodic decrease in U 
describes relaxation of vortices in the pinning wells. With the 
max-field Uo remaining constant over these regions, a reduc¬ 
tion <5-Fpin in pinning force density below F c shows up. 


A. Transient initialization regime 

Our task is to solve the boundary-driven differential 
equation 

-d t U + Dd 2 x U- ^P(C/-£/o) = 0 (10) 

V 

with the diffusion constant D = Bq/Atyt], the max- 
field Uo(X,t) = max t ' <t U(X,t'), and the external drive 
d x ll{0,t ) = —(h ac /B 0 )e~ lult . In a sample of finite thick¬ 
ness d along X vortices cannot move beyond the sample 
center, hence U(d/ 2, t) = 0 provides the second boundary 
condition. We consider a situation where the magnetic 
field Bo changes little over the critical state profile, e.g., 
as it is the case for a fully penetrated sample with thick¬ 
ness d <C L, L = cBo/Anjc the asymptotic extension of 
the critical state profile, B{X) = Bq(L — X)/L. 


Consider first the driven diffusion equation without the 
pinning term. The (asymptotic) solution of this equation 
has the form U{X,t) = (h ac /B 0 )e~^ 1 ~ l ^ x ^ s ' > e~ lujt with S 
the skin depth involving the flux-flow conductivity; the 
displacement U(0,t) at the boundary describes the net 
flux (per unit length along Y) <f>{t) = B o U{0,t) that has 
entered the sample up to time t —obviously flux period¬ 
ically enters and leaves the sample as well known from 
the skin effect. 

The pinning force term — a sp (U — Uq) breaks the sym¬ 
metry between vortices moving in and out of the sample. 
While flux entry proceeds as before, flux exit is inhibited 
as the pinning force pushes the vortices into the sample 
(due to the reduction of the critical force). As a result, 
flux is periodically pumped into the sample until the to¬ 
tal additional flux reaches a value <j> = h ac d/2 and the 
internal magnetic field is shifted from Bq to Bq + h ac . 

The differential equation for the displacement field 
U(X,t) as well as the max-field Uo{X,t) can be found 
by numerical integration of Eq. (10). To this end, we 
bring Eq. (10) into a dimensionless form by measuring 
the time in units of the period r = 2tt/uj and lengths (and 
displacements) in units of the typical diffusion length 
Id = dD t during the period r. The differential equation 
(10) then assumes the form 

-d t U + d 2 x U- < £(U-Uo) = 0, (11) 

with t, X, U, and Uq the new dimensionless quanti¬ 
ties and A c = y/Drj/a sp = \JB^/Aircis-p the strong pin¬ 
ning result for the Campbell length. The boundary con¬ 
dition at the surface X = 0 now reads dxU{0,t) = 
— (h ac /B 0 )e~ l2 ' Kl: , while the vanishing of the displacement 
at the sample center, U{d/2£n,t) = 0, inhibits vortices 
from further penetrating the sample. The relative am¬ 
plitude {h ac /B q) trivially scales the solution, while the 
ratios A c /^d and d/2iD determine the shape of the solu¬ 
tion for a specific setup. The algorithm to solve Eq. (11) 
is based on the forward Euler routine on a discrete mesh 
with spacing AX {At) in the spacial (time) domain. The 
standard scheme 


At P 2 At 

U{X , t + At) = U{X, t) + [U{X - AX, t) - 2U{X, t) + U{X + AX, t)] - [U{X, t) - U 0 {X, <)], (12) 


has to be supplemented with an update of the max-field 

U 0 {X,t +At) =max{U{X,t +At),U 0 {X,t)}. (13) 

This additional rule does not allow for more advanced nu¬ 
merical approaches such as the implicit Crank-Nicolson 
method. The physical constraint to resolve the dynamics 
within the Campbell length (AX) 2 <C (A c /(d) 2 sets a 


lower bound for A A'. The requirement of good conver¬ 
gence of the numerical code imposes a constraint on the 
mesh size At via At/{AX ) 2 1 [these two constraints 

also guarantee that ^At/A 2 1, see Eq. (12)]. A sec¬ 
ond constraint on the time mesh At <C 1, guaranteeing 
the smoothness of the external drive, is automatically 
satisfied for the parameter range considered here. The re- 
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suits shown in Figs. 3-5 assume X c /£n = l/v / 207r m 1/8 
and d/2(-D = 1. In Fig. 6 different system sizes d/2 £d = 
1, 3, 5, 7, and 9 are shown. 

The displacement fields U(X, t ) and Uo(X, t) are shown 
in Fig. 3 for several positions X, once at the boundary 
X = 0 where U is proportional to the penetrated flux <f> 
(per unit length along Y), as well as a few penetration 
depths A c into the sample at X = 2X C and X = 4A C . 
Due to the action of the pinning force term, the dis¬ 
placement U is no longer periodic: the increase in U 
describing flux entry is interrupted by short excursions 
with decreasing U as vortices relax in their pinning po¬ 
tential. The max-field Uq ignores these depressions and 
the finite amplitude Uq — U generates the decrease in 
pinning force SF pin . The latter is shown in Fig. 4: Re¬ 
gions of vanishing SF pul where additional flux penetrates 
deeper into the sample are interrupted by segments of 
finite force SF p i n > 0 where vortices relax back in their 
pinning potentials. With increasing time, the additional 
flux shifts the internal field upwards to B 0 + h ac . For 
long times, vortices evolve reversibly within their pin¬ 
ning wells, with no further flux entering the sample. As 
a result, the change in pinning force density SF p - ln os¬ 
cillates back and forth with maximal amplitude and no 
intermediate regions with SF p - ln = 0 show up. 

The physically most relevant and transparent result 
is the behavior of the magnetic induction B within the 
sample. In Fig. 5 we show the evolution of SB(X , t) with 
time t over several penetration depths into the sample 
along X. Additional flux is periodically pumped into the 
sample, lifting the dc value of B with increasing time t 
(change from green to blue color). For better illustration 
we have chosen a geometry and parameters such that 




FIG. 5. Change in magnetic induction 8B(X,t) as a function 
of time t and distance X into the sample. Every oscillation 
pumps flux into the sample, pushing the critical state at Bo 
(SB = 0, green) towards one at Bo + h ac (SB = h ac , blue). 
The 8 cycles shown here push almost all the flux into the 
sample that is required for its shift by h ac - 

the 8 cycles shown in the figure nearly suffice to reach 
the asymptotic value B 0 + h ac , see also the lowest curve 
in Fig. 6 for d/2W = 1. 

In the following, we build and analyze a simple model 
for the flux entry into the sample; this will allow us to 
estimate the number of cycles needed to reach the asymp¬ 
totic state. We start with the first ac cycle which pumps 
a flux (per unit transverse length) 4>o = f h ac i d into the 
sample, where £d = VDt is the diffusion length during 
one cycle, r = 2i r/cu, and / < 1 is a numerical factor ac¬ 
counting for the magnitude of penetrated flux per pulse 
(note that flux is pumped into the sample only during 
the first half-cycle and the precise magnitude depends 
on the detailed shape of the driving signal). This flux 
pulse diffuses with time t and spreads over the distance 
\J~Dt. At time t, the flux pulse then reduces the ac- 
magnetic field amplitude at the sample boundary X = 0 
by ho(t) = cj>o/y/Dt. Hence the next pulse <j>i entering 
the sample is reduced, (f>i = f [h ac — Tio(t)\£d- Iterating 
this process, the i-th flux pulse entering the sample is 
given by 

4>i = fhi £ D , (14) 

with the iteratively defined amplitude 


FIG. 4. Change in pinning force density SF P i n versus dis¬ 
placement field U at the sample surface X = 0. The time 
t is a curve parameter that evolves as indicated by the red 
line. Starting from U = 0 and increasing U , flux moves into 
the sample, the pinning force density is at its maximum value 
F/in = — F c , and 5F pln = 0. Deviations of F p j n from its crit¬ 
ical value appear whenever U decreases and vortices relax in 
their pinning potential. At long times, the Bean profile has 
been raised to the value Bo + h ac and vortices oscillate re¬ 
versibly (see double-headed red arrow) within their potential 
wells without additional flux entrance and vanishing regions 
with SFpi n = 0. The max-field Uo(0,t ) reaches its asymptotic 
value Uq(X = 0), see Eq. (24). 


and 


hi — hac 


T, h m (ir ) 

771<i 


(15) 


h m (t) = <\>ml\/D(t - tot). (16) 


Combining Eqs. (14), (15), and (16) we then have to solve 
the self-consistency equation for hi, 


hi — h ac f ^ ' 


\/i — : 


(17) 
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Going over to continuous variables, we obtain the integral 
equation 

h{t) = h ac - -£= [ dt'-^==, (18) 

V T Jto vt - t' 

with the starting time to to be determined together with 
h(t). Inserting the Ansatz h(t) = jd h ac /y/i, we can carry 
out the integral and arrive at the self-consistency condi¬ 
tion 


P 


hac 

Vt 



vi 


arcsin (y/v/i) 



(19) 


Choosing /3 = yfrjf V (such that the upper boundary 
t provides a term (2/n) arcsin 1 = 1) and setting to = 
t/2/ 2 , we obtain a consistent solution for large times 1 ' 
t to in the form 


(20 > 

Finally, adding up the flux pulses (j>i, we find that the 
flux penetrates into the sample diffusively as 

4>{t) = -hadoWi ( 21 ) 

7T 

independent on /. The flux needed to push the internal 
critical state from an initial field Bod/2 to a final field 
(Bo + h ac )d/2 then involves the time 


_ Vd^ 

= is tl T 


( 22 ) 


or n = (ir 2 /16) d 2 /£ 2 D cycles. In Fig. 6 we compare the 
result of our analytic model calculation with our numer¬ 
ical data. For the latter, we have computed the response 
of the vortex system for different system sizes with the 
thickness d/2 of the critical state profile taking the values 
d/2 £d = 1, 3, 5, 7, and 9 in units of the diffusion length 
£d- Despite the simplicity of the model, the analytic re¬ 
sult in Eq. (21) is in good qualitative agreement with 
the numerical results; for a quantitative agreement, the 
model result has to be multiplied with a factor « 1.25. 

Our comparison with recent experiments 9 refers to a 
SrPd 2 Ge 2 single crystal, a material that is isostructural 
to the Fe- and Ni-pnictides, of typical size d ~ 0.5 mm in 
a field H ~ 0.1 T and subject to an ac field h ac = 2 pT 
at frequencies v = w/27t « 17 MHz. With a normal 
state resistivity 18 p n ~ 20 /iflcm and the upper crit¬ 
ical field 9 H c2 « 0.5 T, we obtain a diffusion length 
£d ~ 40 /mi, hence d/2£jj ~ 6. The pinning strength is 
quantified through the Campbell penetration depth mea¬ 
suring A c ~ 10 /mi, hence A c /£d ~ 1/4. The elementary 
displacement 17(0, r) ~ £d(h ac /H) is of the order of 8 A 
at the beginning of the initial diffusion dominated regime 
and smaller, £7(0, r) ~ X c (h ac /H) ss 2 A in the asymp¬ 
totic regime, far below the coherence length £ ss 25 A. For 
such a sample and with the above experimental parame¬ 
ters we find that the asymptotic state is reached within 



FIG. 6. Flux pumping <j>(t) = BqU (0, t) as a function of time. 
We compare our numerical results (oscillating curves) for dif¬ 
ferent sample sizes d with our analytical model (dotted line), 
see Eq. (21). A better quantitative agreement with the nu¬ 
merical data is obtained when scaling the analytic result by 
1.25 (dashed line). The horizontal lines indicate for each sys¬ 
tem the value of the penetrated flux h ac d/2 after termination 
of the initialization process. The inset shows an expanded 
view at small times of both U and Uo at the surface. 


about one hundred cycles or ij n j t ~ 10 /ts, i.e., the tran¬ 
sient initialization regime is usually not interfering with 
the measurement of the Campbell length A c . Each of 
these cycles typically pumps a fraction of a vortex into 
the sample, as (<£od)/$o ~ 1/5. To complete the dis¬ 
cussion, we briefly describe the result for the asymptotic 
regime. 


B. Asymptotic periodic regime 

Once the asymptotic time domain has been reached, 
vortices exhibit the typical oscillatory behavior near the 
surface, but with respect to the new critical state that 
has been shifted upward by h ac . Eq. (3) can be solved 
exactly for a half-infinite space and its solution 

U(X, t ) = U 0 (X) - \ c (h ac /Bo)e~ x / Xc [1 - e~ iuit ) (23) 

with the asymptotic max-field, 

U 0 (X) = (cj>-h ac X)/Bo, (24) 

and the penetrated dc flux <j> provides a good approxima¬ 
tion for a finite-size sample with d/X c 1. In Eq. (23), 
the first term Uo describes an overall shift of the crit¬ 
ical state profile, SBq = —Bod x Uo = h ac , and hence 
4> = hacd/2. The second term in Eq. (23) captures 
the oscillatory and decaying part of the magnetic field 
—h ac e~ x l Xc [ 1 — e - ’"*]. Alternatively, the asymptotic 
change of the magnetic field profile can be decomposed 
into an ac component SB ac (X,t ) = h ac e~ x / Xc e~ lult and 
an inhomogeneous dc component 6Bd c (X) = h ac [l — 
e~ x ' Xc ], The conversion of the ac drive into a dc sig¬ 
nal is a unique feature of the critical state and may be 
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observed in an experiment. The Campbell penetration 
depth A c is determined by a sp , Eq. (9), involving the 
jump in force A/ pin , 

A C = ^-(25) 

47ra sp ripCL^t ; z k 

with k the Labusch parameter, k ps / p ao/eoC- Here, f p 
denotes the pinning force and £q = ($o/47rA L ) 1 2 3 4 5 6 is the 
vortex line energy. The small parameter n p aQ^ 2 K <C 1 
defines the three-dimensional strong pinning regime 14 . 

The jump in force Af pin in general depends on the 
state preparation of the vortex system 10 . Hence field- 
cooled and zero-field cooled vortex states may exhibit 
different penetration depths and hysteretic effects may 
show up. This better microscopic understanding of the 
Campbell penetration depth may then allow to gain more 
detailed information on the pinscape. 

C. Non-linear ac-response 

The above discussion has focused on Campbell pene¬ 
tration, i.e., the linear response regime at small drive h ac - 
As mentioned in the introduction, this has to be distin¬ 
guished from the Bean penetration at large amplitudes 
h ac • The two regimes are separated by the conditions 
Sj ~ jc or h ac ~ JcAc/c, where the vortex displacement 
U ~ (h ac /B q)X g matches the pinning length x+ ~ 
and k > 1 is the strong pinning parameter 12 . The Bean 
penetration at large amplitudes exhibits similar behavior 
on a qualitative level. E.g., using Bean’s original (qua¬ 
sistatic) approach, the dc shift of the critical state profile 
by h ac occurs during the first (quarter of the) ac cycle; 
making use of a more sophisticated model with a specific 
I-V characteristics would result in a penetration involv¬ 
ing flux-pulses similar to those found above. At large 
times, after the critical state has been shifted upwards 
by the ac amplitude h ac , the Bean penetration reaches 
an oscillatory regime where the ac drive generates a sim¬ 
ple periodic dynamics on the penetration scale, see Fig. 1 
(right panel). On the other hand, the Bean penetration 
differs from the Campbell penetration on a quantitative 
level: the penetration depth ~ ch ac /jc within the 
Bean model scales with h ac , resulting in a third harmonic 


ac magnetic response 19 . The large displacements (com¬ 
pared to the pinning length x+) in the Bean regime do 
not resolve the internal structure of the pinning centers 
and hence the ac magnetic response is independent of the 
state preparation. Finally, the Bean penetration exhibits 
generic hysteretic properties. 


IV. SUMMARY AND CONCLUSION 

Using the results of strong pinning theory, we have 
studied the complete ac magnetic penetration dynamics 
in the Campbell regime for a superconducting sample 
prepared in a critical state. The field penetration pro¬ 
ceeds in two phases, a short (on the time scale of the ex¬ 
periment) transient initialization regime where flux en¬ 
ters the sample in a sequence of pulses that shifts the 
overall critical state by the ac amplitude h ac , followed by 
the standard ac dynamical regime where vortices merely 
move back and forth in their pinning potentials with no 
additional net flux entering the sample. 

The apparent divergence of the Campbell penetration 
depth A 0 in the phenomenological description 8 leaves its 
trace in the pulsed diffusive flux penetration through¬ 
out the sample during the initialization regime. After 
long times, the Campbell length A c is a (regular) linear- 
response parameter containing valuable information on 
the pinscape that can be extracted from its dependence 
on the state preparation 10 . Despite qualitative similar¬ 
ities between the Campbell (at low ac fields) and the 
Bean penetration (at higher ac fields) the latter differs 
by i) its field-dependent penetration depth £ B ~ ch ac /jc , 
ii) the generation of a third harmonics ac signal and iii) 
the insensitivity to the pinscape structure. Finally, the 
rectified dc signal 6Bd c (X) induced by the ac (Campbell 
or Bean) drive is a signature that can be probed in future 
experimental investigations. 
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